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1. Introduction 


Let X],..., X n be n > 1 independent copies of a non-negative random variable 
(rv) X, defined over some probability space (12, M, P), with absolutely continuous 
cumulative distribution function (cdf) F. An independent sequence of independent 
rv’s Y 1 ,...,Y n , with absolutely continuous cdf G, censor them to the right, so that 
at each stage j we can only observe Zj := min(Xj, Yj) and 8j := 1 {Xj < Yj} , with 
1 {•} denoting the indicator function. The rv Sj indicates whether or not there has 
been censorship. Throughout the paper, we use the notation S(x) := S(oo) —S(x), 
for any S. If H denotes the cdf of the observed Z's, then, by the independence of 
Xi and Y l , we have H (z) = F (z) G ( z ). In our work, we assume that both F and G 
are heavy-tailed, this means that there exist to constants 71 > 0 and 72 > 0 , called 
tail indices, such that 


lim = *-!/■» and lim 3A = *-!/■», 

z-KX) F(z) z^oo G(z) 


( 1 . 1 ) 


for any x > 0. Consequently, H is heavy-tailed too, with tail index 7 := 7172/(71 + 
72 ). The class of heavy-tailed distribution takes a significant role in extreme value 
theory. It includes distributions such as Pareto, Burr, Frechet, a—stable (0 < a < 2) 
and log-gamma, known to be appropriate models for fitting large insurance claims, 
log-returns, large fluctuations of prices, etc. (see, e.g., Resnick, 2007). Examples 
of censored data with apparent heavy tails can be found in Gomes and Neves (2011). 
The nonparametric maximum likelihood estimator of F is given by Kaplan and Meier 
(1958) as the product limit estimator 


F n (x) 


1- n 


n-j 
n — j + 1 


%:n] 


for x < Z n:n 
for x > Z n:n , 


where Z 1:n < ... < Z n:n denote the order statistics pertaining to the sample (Z 1 ,..., Z n ) 
with the corresponding concomitants <5[i :n ], •••, d[ n:n ] satisfying 5[j :n ] = Si if Zj :n = Zi. 
This estimator, known as Kaplan-Meier estimator of F, may be expressed as follows 


F n (x) := y; W i>n 1 {Z i:n < x} , (1.2) 

i=2 


where, for 2 < i < n, 


Wi n := 


n — i + 1 


i— 1 

n 

3 = 1 


71-J 


%:n] 


n — j + 1 
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(see, e.g., Reiss and Thomas, 2007, page 162). The aim of this paper is to propose 
an asymptotically normal estimator for the mean of A", 


/j := E[A] = / F(x)dx , 


whose existence requires that 71 < 1. The sample mean for censored data is obtained 
by substituting, in the previous equation, the cclf F by its estimator F n to have 


p : = 




r *-1 

TT 

n — i + 1 1 

i= 2 3=1 


n-j 
n — j + 1 


S[j:n] 


The asymptotic normality of p n is established by Stute (1995), under the assump¬ 
tions that the integrals 


h := l and /, := / X ' dF(x), 

be finite, where T 0 (x) := exp {JJ dH^(z)/H(z)J with (v) := P (Z < v, 5 = j ), 
j = 0,1. In the sequel, the latter functions will play a prominent role. How¬ 
ever, when we deal with heavy-tailed distributions, the quantities I\ and I 2 may 
be infinite. Indeed, suppose that both F and G are Pareto distributions, that is 
F(x) = x~ lhl and G(x) = x~ 1 ^ 2 , for x > 1. This obviously gives H(x) = x -1 / 7 , 
H^){x) = 7(1 — x - 1 /, 7 )/ 72 , H^ l \x) = 7(1 — x~ l ^)/^\ and T 0 (x) = x 1 ^ 2 . Whenever 
( 71 , 72 ) are such that 71 > 72 / (1 + 272 ), we readily check that I\ — I 2 — 00 . I 11 
other words, the range 


K ■= < 7 i >72 > 0 : 


72 


< 71 < ! 


1 + 272 

is not covered by the central limit theorem established by Stute (1995), and thus, 
another approach to handle this situation is needed. This problem was already 
addressed by Peng (2001) for sets of complete data from heavy-tailed distributions 
with tail indices lying between 1/2 and 1. Note that in the non censoring case, we 
have 71 = 7 meaning that 72 = 00 , consequently 1Z reduces to Peng’s range. The 
consideration of the range 1Z is motivated and supported from a practical point of 
view as well. Indeed, as an example Einmahl et al. (2008) analyzed the Australian 
AIDS survival dataset and found that 71 = 0.14 and p = 0.28 leading to 72 = 0.05. It 
is easily checked that these index values belong to TZ and therefore Stute’s result does 
not apply in this situation. To define our new estimator, we introduce an integer 
sequence k = k n , representing a fraction of extreme order statistics, satisfying 


1 < k < n, k —> 00 and k/n —* 0 as n —» 00 , 


(1.3) 
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and we set h = h n := F _1 (l — k/n), where K := inf { x : K(x) > y} ,0 < y < 1, 
denotes the quantile function of a cdf K. We start by decomposing p as the sum of 
two terms as follows: 

ph poo 

p = / F{x)dx + / F(x)dx =: p\ + p 2 , 

Jo Jh 

then we estimate each term separately. Integrating the first integral by parts and 
changing variables in the second respectively yield 

— r h — f°° F (hr) 

Pi = hF(h) + / xdF(x) and /i 2 = hF (h) / —- -dx. 

Jo J1 F (h) 

By replacing h and F(x) by Z n _k-.n and F n (x) respectively and using formula (1.2) , 


we get 


hi := 


n—k 

n 

3 =1 


n-j 
n — j + 1 


n—k 


J n—k 


:n + Yl 


2—1 

n 


n — i + 1 In - 7 + I 

i =2 i=l V J 


n ~ J 




7 


(1.4) 


as an estimator to /j j . Regarding /x 2 , we apply the well-known Karamata theorem 
(see, for instance, de Haan and Ferreira, 2006, page 363), to write 

7i 


h 2 


7i 


-hF (h), as n —> 00 , 0 < 71 < 1 . 


The quantities h and F (h) are, as above, naturally estimated by Z n _ k . n and 


F (Z n _ k:n ) — JJ 


n—k 


n-j 


L i=i \n — j + 1 

respectively. Now, it is clear that to derive an estimator to yU 2 , one needs to estimate 
the tail index 71 . The general existing method, which first appeared in Beirlant et al. 
(2007) and then developed in Einmahl et al. (2008), is to consider any consistent 
estimator of the extremal index 7 based on the Z-sample and divide it by the 
proportion of non-censored observations in the tail. For instance, Einmahl et al. 
(2008) adapted Hill’s estimator to introduce an estimator 7 ^'^ := 7 H /p to the tail 
index 7 x under random right censorship, where 

-i ^ ry 1 ^ 


'■= 7 l0g ^7 l+l n and P '= T S ^H-i+l:n], 

k *7-1 Z n — k:n k *—> 


1 —it — rb. it • -1 

2=1 2=1 

with k = k n satisfying (1.3) , are the classical Hill estimator and the proportion of up¬ 
per non-censored observations respectively. It is proved in Brahimi et al. (2015) that 
p consistently estimates p := 72 / (71 + 72 ), therefore 'y[ H ' c> consistently estimates 
71 = 7 Jp. The authors of Brahimi et al. (2015) provide a Gaussian approximation 
leading to the asymptotic normality of 7 ^’^ by adopting a different approach from 
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that of Einmahl et al. (2008), who also showed that is asymptotically normal. 
Consequently, we obtain 


H,c ) n—k 

- n ry TT l n — J 

^2 ■- Z ™-fc:n II 


n — j + 1 


+ :™] 


, for 7 ^’ c) < 1 , 


(1.5) 


i - 7** 

as an estimator to /a 2 . Finally, with (1.4) and (1.5) , we construct our estimator fi of 
the mean /i : 


n—k r i —1 / . \ 

<2 . = TT ( n ~ X 

4-^ n — i + lJ-J- In — 7 + 1 
7=2 j=l V J 


n—k / \ S 

n — J 


n 

1=1 


n — j + 1 


[jin] 7 

ZJ r. 


n—k:n 


1 - 7f’ c) ' 


The rest of the paper is organized as follows. In Section 2, we state our main result 
which we prove in Section 4. Section 3 is devoted to a simulation study in which we 
investigate the finite sample behavior of the newly proposed estimator J2. Finally, 
some results, that are instrumental to our needs, are gathered in the Appendix. 

2. Main results 

Our main result, established in the following theorem, consists in the asymptotic 
normality of the newly introduced estimator j2. We notice that the asymptotic nor¬ 
mality of extreme value theory based estimators is achieved in the second-order 
framework (see de Haan and Stadtmiiller, 1996). Thus, it seems quite natural to 
suppose that cdf’s F and G satisfy the well-known second-order condition of regular 
variation. That is, we assume that there exist constants Tj < 0 and functions Aj, 
j = 1,2 tending to zero, not changing sign near infinity and having regularly varying 
absolute values with indices Tj, such that for any x > 0 


lim 

£—>■00 


F(tx)/F(t)—x 1//71 


Ai(f) 


— j.. I/u 


X 


.71 /71 


_ ^.-1/72 


[im G(te)/G(t) - x ■- =x _ lh , 

£—>00 


7i r i 

272/72 _ 1 


( 2 . 6 ) 


A 2 (t) 72 r 2 

Theorem 2.1. Assume that the second-order conditions of regular variation (2.6) 
hold with 72 / (1 + 272 ) < 71 < 1. Let k = k n be an integer sequence satisfying, in 
addition to (1.3), lim ra _>. 0 0 \/kAi(h) < 00 and \fkhF (h) —> 00 . Then there exist 
finite constants m and a 2 > 0 such that 


Vk(fl- n) d 


^n—k-.nFn^^n—k-.nf 

Remark 2.1. We have 

Ai 


m := 


4 Af( m, cr 2 ) , as n —> 00 . 


Ai 


(1 -pTi) (1 - 7i) 2 (71 + n - 1) (! - 7i)’ 
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with Ai := lim^oo \/kAi(h), whereas the computations of the asymptotic variance 
a 2 are very tedious and result in an expression that is too complicated. However, 
the lack of a closed form for a 2 could be overcome in applications, as both parame¬ 
ters are usually estimated by the respective sample mean and variance obtained by 
bootstrapping fi. 


3. Simulation study 

We carry out a simulation study to illustrate the performance of our estimator, 
through two sets of censored and censoring data, both drawn, in the first part, from 
Frechet model 


F (x) = exp {—x 71 } , G(x) — exp {—x 72 } , x > 0, 

and, in the second part, from Burr model 

F (x) = 1 - (1 + x 1/v y vhl , G (x) = 1 - (1 + *^p /72 , X > 0, 

where 77,71,72 > 0. We fix 77 = 1/4 and choose the values 0.3, 0.4 and 0.5 for 7!. 
For the proportion of the really observed extreme values, we take p = 0.40, 0.50, 
0.60 and 0.70. For each couple (71, p ), we solve the equation p = 72/(71 +72) to get 
the pertaining 72-value. We vary the common size n of both samples (Ad,..., X n ) 
and (Yi, ...,Y n ), then for each size, we generate 1000 independent replicates. Our 
overall results are taken as the empirical means of the results obtained through the 
1000 repetitions. To determine the optimal number (that we denote by k*) of upper 
order statistics used in the computation of 7^ , we apply the algorithm given 
in page 137 of Reiss and Thomas (2007). The performance of the newly defined 
estimator fi is evaluated in terms of absolute bias (abs bias), mean squared error 
(mse) and confidence interval (conf int) accuracy via length and coverage probability 
(cov prob). 

The results, summarized in Tables 3.1, 3.2 and 3.3 for Frechet model and Table 3.4, 
3.5 and 3.6 for Burr distribution, show that the same conclusions might be drawn in 
both cases. As expected, the sample size influences the estimation in the sense that 
the larger n gets, the better the estimation is. On the other hand, it is clear that 
the estimation accuracy increases when the censoring percentage decreases, which 
seems logical. Moreover, the estimator performs best for the smaller value of the 
tail index, as we can see from Tables 3.1 and 3.4. Finally, many simulations realized 
with extreme value indices larger than 0.5, but whose results are not reported here, 
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yi — 0.3 —y ii — 1.298 


p = 0.40 


n 

P 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.247 

0.051 

0.021 

1.043 - 

1.450 

0.88 

0.407 

1000 

1.244 

0.054 

0.020 

1.099 - 

1.389 

0.88 

0.291 

1500 

1.233 

0.065 

0.005 

1.119 - 

1.346 

0.80 

0.227 

2000 

1.231 

0.067 

0.005 

1.135 - 

1.328 

0.74 

0.193 




P - 

= 0.50 




500 

1.248 

0.050 

0.008 

1.049 - 

1.447 

0.96 

0.399 

1000 

1.247 

0.051 

0.004 

1.107- 

1.387 

0.90 

0.280 

1500 

1.250 

0.048 

0.003 

1.134 - 

1.365 

0.90 

0.231 

2000 

1.248 

0.050 

0.003 

1.146 - 

1.350 

0.86 

0.204 




P - 

= 0.60 




500 

1.254 

0.044 

0.009 

1.050 - 

1.458 

0.90 

0.408 

1000 

1.257 

0.041 

0.003 

1.119 - 

1.395 

0.94 

0.275 

1500 

1.266 

0.032 

0.002 

1.153 - 

1.379 

0.96 

0.226 

2000 

1.264 

0.034 

0.002 

1.164 - 

1.364 

0.92 

0.200 




P - 

= 0.70 




500 

1.265 

0.033 

0.003 

1.069 - 

1.460 

0.97 

0.391 

1000 

1.269 

0.029 

0.002 

1.123 - 

1.415 

0.96 

0.291 

1500 

1.279 

0.019 

0.001 

1.162 - 

1.395 

0.98 

0.233 

2000 

1.278 

0.020 

0.001 

1.178 - 

1.377 

0.96 

0.199 


Table 3.1. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Frchet model with shape parameter 0.3 


show that the estimator behaves poorly especially when the censorship proportion 
is high. 



7 i = 0.4 —y n =1.489 


p = 0.40 


n 

p 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.370 

0.120 

0.074 

1.147- 

1.593 

0.71 

0.446 

1000 

1.377 

0.112 

0.048 

1.217- 

1.536 

0.57 

0.319 

1500 

1.367 

0.122 

0.019 

1.241 - 

1.493 

0.48 

0.252 

2000 

1.363 

0.126 

0.018 

1.256 - 

1.470 

0.36 

0.214 




P z 

= 0.50 




500 

1.396 

0.093 

0.027 

1.169 - 

1.624 

0.81 

0.455 

1000 

1.394 

0.095 

0.018 

1.237- 

1.551 

0.66 

0.313 

1500 

1.392 

0.097 

0.012 

1.264 - 

1.521 

0.65 

0.257 

2000 

1.389 

0.101 

0.012 

1.275 - 

1.502 

0.55 

0.227 




P z 

= 0.60 




500 

1.407 

0.082 

0.013 

1.189 - 

1.625 

0.89 

0.436 

1000 

1.405 

0.084 

0.010 

1.251 - 

1.559 

0.77 

0.308 

1500 

1.419 

0.070 

0.007 

1.292 - 

1.546 

0.84 

0.254 

2000 

1.418 

0.071 

0.007 

1.308 - 

1.529 

0.71 

0.222 




P z 

= 0.70 




500 

1.420 

0.069 

0.010 

1.199 - 

1.641 

0.92 

0.442 

1000 

1.433 

0.056 

0.006 

1.273 - 

1.593 

0.86 

0.320 

1500 

1.443 

0.046 

0.004 

1.312 - 

1.575 

0.90 

0.263 

2000 

1.442 

0.047 

0.004 

1.329 - 

1.554 

0.89 

0.226 


Table 3.2. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Frchet model with shape parameter 0.4 


4. Proofs 

We begin by a brief introduction on some uniform empirical processes under random 
censoring. The empirical counterparts of (j — o, 1) are defined, for v > 0, by 

1 n 

H n ] (v) :=-Vl {Zi< v, Si = j} , j = 0 , 1 . 
n 

i =1 
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7 i = 0.5 -+n = 1.772 


p = 0.40 


n 

P 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.566 

0.206 

0.398 

1.262 - 

1.870 

0.52 

0.608 

1000 

1.550 

0.223 

0.176 

1.372 - 

1.727 

0.28 

0.355 

1500 

1.559 

0.214 

0.064 

1.415 - 

1.703 

0.20 

0.289 

2000 

1.549 

0.224 

0.061 

1.426 - 

1.671 

0.13 

0.245 




P z 

= 0.50 




500 

1.577 

0.195 

0.180 

1.309 - 

1.846 

0.53 

0.537 

1000 

1.573 

0.199 

0.139 

1.386 - 

1.761 

0.37 

0.375 

1500 

1.578 

0.195 

0.051 

1.430 - 

1.725 

0.20 

0.294 

2000 

1.576 

0.196 

0.044 

1.447 - 

1.706 

0.22 

0.259 




P z 

= 0.60 




500 

1.626 

0.147 

0.128 

1.362 - 

1.889 

0.65 

0.527 

1000 

1.617 

0.155 

0.034 

1.430 - 

1.805 

0.56 

0.375 

1500 

1.606 

0.166 

0.033 

1.465 - 

1.747 

0.34 

0.282 

2000 

1.622 

0.150 

0.029 

1.494 - 

1.751 

0.34 

0.258 




P z 

= 0.70 




500 

1.632 

0.141 

0.046 

1.375 - 

1.888 

0.72 

0.513 

1000 

1.646 

0.126 

0.024 

1.459 - 

1.833 

0.70 

0.370 

1500 

1.668 

0.104 

0.017 

1.516 - 

1.821 

0.68 

0.305 

2000 

1.666 

0.107 

0.016 

1.535 - 

1.797 

0.57 

0.262 


Table 3.3. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Frchet model with shape parameter 0.5 


In the sequel, we will use the following two empirical processes 
Vn (Hn\v) ~ H U) (v ;)) , j = 0,1; v > 0, 

which may be represented, almost surely, by a uniform empirical process. Indeed, 
let us define, for each i = 1, ...,n, the following rv 


Ui := 5iHW(Zi) + (1 - 5i)(e + H^(Zi)). 
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yi — 0.3 —y fi — 1.228 


p = 0.40 


n 

V> 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.186 

0.042 

0.077 

0.972 - 

1.399 

0.90 

0.428 

1000 

1.179 

0.049 

0.019 

1.038 - 

1.32 

0.80 

0.282 

1500 

1.163 

0.064 

0.005 

1.053 - 

1.273 

0.80 

0.220 

2000 

1.164 

0.063 

0.005 

1.068 - 

1.261 

0.72 

0.193 




P 

= 0.50 




500 

1.186 

0.042 

0.009 

0.991 - 

1.380 

0.94 

0.388 

1000 

1.173 

0.054 

0.004 

1.039 - 

1.308 

0.93 

0.269 

1500 

1.068 

0.047 

0.003 

1.180 - 

1.292 

0.88 

0.224 

2000 

1.181 

0.046 

0.003 

1.086 - 

1.276 

0.86 

0.190 




P 

= 0.60 




500 

1.184 

0.043 

0.004 

0.997 - 

1.371 

0.95 

0.374 

1000 

1.192 

0.036 

0.002 

1.058 - 

1.326 

0.96 

0.268 

1500 

1.196 

0.031 

0.002 

1.088 - 

1.305 

0.96 

0.217 

2000 

1.194 

0.034 

0.002 

1.099 - 

1.288 

0.92 

0.190 




P 

= 0.70 




500 

1.198 

0.029 

0.003 

1.012 - 

1.384 

0.97 

0.373 

1000 

1.200 

0.028 

0.001 

1.066 - 

1.334 

0.98 

0.269 

1500 

1.208 

0.020 

0.001 

1.098 - 

1.317 

0.98 

0.219 

2000 

1.207 

0.021 

0.001 

1.113 - 

1.301 

0.98 

0.188 


Table 3.4. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Burr model with shape parameter 0.3 


From Einmahl and Koning (1992), the rv’s U \,..., U n are iid (0, l)-uniform. The em¬ 
pirical cdf and the uniform empirical process based upon U\, ...,U n are respectively 
denoted by 

1 n 

U n (s) : =- 1 {Ui < s} and a n (s) : = V^(U„(s) - s), 0 < s < 1. 

n 

i —1 

Deheuvels and Einmahl (1996) state that almost surely 

H®(v) = U n (H(°\v) + 9)- U n (0), for 0 < H^\v) <1-9, 
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7i = 0.4 —» ji = 1.498 


p = 0.40 


n 

h 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.426 

0.071 

0.093 

1.193 - 

1.660 

0.76 

0.466 

1000 

1.388 

0.110 

0.033 

1.224 - 

1.551 

0.58 

0.327 

1500 

1.374 

0.124 

0.020 

1.248 - 

1.499 

0.44 

0.252 

2000 

1.374 

0.123 

0.019 

1.268 - 

1.480 

0.29 

0.212 




P z 

= 0.50 




500 

1.402 

0.096 

0.047 

1.176 - 

1.627 

0.80 

0.451 

1000 

1.389 

0.109 

0.017 

1.231 - 

1.546 

0.64 

0.316 

1500 

1.401 

0.097 

0.012 

1.272 - 

1.530 

0.66 

0.258 

2000 

1.402 

0.096 

0.011 

1.292 - 

1.511 

0.53 

0.219 




P z 

= 0.60 




500 

1.422 

0.076 

0.043 

1.186 - 

1.657 

0.85 

0.471 

1000 

1.421 

0.077 

0.009 

1.261 - 

1.581 

0.86 

0.320 

1500 

1.429 

0.069 

0.007 

1.302 - 

1.556 

0.80 

0.254 

2000 

1.427 

0.071 

0.006 

1.316 - 

1.538 

0.76 

0.223 




P z 

= 0.70 




500 

1.436 

0.061 

0.009 

1.214 - 

1.658 

0.94 

0.444 

1000 

1.441 

0.057 

0.006 

1.285 - 

1.597 

0.92 

0.312 

1500 

1.451 

0.047 

0.004 

1.322 - 

1.580 

0.91 

0.259 

2000 

1.449 

0.049 

0.004 

1.340 - 

1.558 

0.88 

0.218 


Table 3.5. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Burr model with shape parameter 0.4 


and 


H^\v) = U n (i/ ( 1 ) (iO), for 0 < H {1 \v) < 9. 

It is easy to verify that almost surely 

A, (i>) := y/n - H W ( t.)) = a n (D)-a n (« - if ( 1 ) (i>)) , for 0 < H m (v) < 6, 

(4.7) 
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7l = 0.5 -A [i = 1.854 


p = 0.40 


n 

V> 

abs bias 

mse 

conf int 

cov prob 

length 

500 

1.654 

0.200 

0.760 

1.330 - 

1.978 

0.50 

0.649 

1000 

1.648 

0.206 

0.114 

1.460 - 

1.836 

0.26 

0.375 

1500 

1.630 

0.224 

0.098 

1.478 - 

1.782 

0.14 

0.304 

2000 

1.621 

0.233 

0.090 

1.491 - 

1.752 

0.14 

0.260 




P 

= 0.50 




500 

1.603 

0.252 

0.554 

1.253 - 

1.952 

0.67 

0.700 

1000 

1.658 

0.196 

0.090 

1.470 - 

1.847 

0.34 

0.378 

1500 

1.653 

0.202 

0.049 

1.501 - 

1.804 

0.25 

0.303 

2000 

1.656 

0.198 

0.045 

1.530 - 

1.782 

0.22 

0.252 




P 

= 0.60 




500 

1.688 

0.166 

0.066 

1.417- 

1.959 

0.67 

0.542 

1000 

1.693 

0.161 

0.036 

1.508 - 

1.879 

0.54 

0.371 

1500 

1.695 

0.159 

0.031 

1.544 - 

1.846 

0.39 

0.301 

2000 

1.705 

0.149 

0.027 

1.576 - 

1.834 

0.34 

0.258 




P 

= 0.70 




500 

1.737 

0.117 

0.060 

1.462 - 

2.012 

0.77 

0.550 

1000 

1.737 

0.117 

0.036 

1.547 - 

1.927 

0.74 

0.380 

1500 

1.749 

0.105 

0.016 

1.593 - 

1.904 

0.70 

0.311 

2000 

1.753 

0.101 

0.014 

1.621 - 

1.885 

0.60 

0.264 


Table 3.6. Absolute bias, mean squared error and 95%-confidence 
interval accuracy of the mean estimator based on 1000 right-censored 
samples from Burr model with shape parameter 0.5 


and 

/3 n ( v ) := \fn (h^\v) — = — a n ^1 — , for 0 < H^°\v) <1 — 9. 

(4.8) 

Our methodology strongly relies on the well-known Gaussian approximation given 
in Corollary 2.1 by Csorgo et al. (1986). It says that: on the probability space 
(G, A, P), there exists a sequence of Brownian bridges { B n (s ); 0 < s < 1} such that 
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for every 0 < ( < 1/4, 


sup 

l/n<s<l 


n c \a n (l - s) - B n ( 1 - s) 
,sV2-C 


— Op(l)- 


(4.9) 


For the increments a n {9) — a n (6 — s), we will need an approximation of the same 
type as (4.9). Following similar arguments, mutatis mutandis, as those used in the 
proofs of assertions (2.2) of Theorem 2.1 and (2.8) of Theorem 2.2 in Csorgo et al. 
(1986), we may show that, for every 0 < 6 < 1 and 0 < ( < 1/4, we have 

n c \{a n (6) - a n {6 - s)} - {B n (9) - B n {9 - s)}[ = ^ ^ 


sup 

1 /n<s<6 


T/2-C 


4.1. Proof of Theorem 2.1. Observe that Ji — // = (// L — fix) + (/j 2 — H 2 ), where 


Hi ~ Hi = 


F n (x)dx — / F(x)dx , 


and 


n—k 


H2 H2 ~ I] ( 1 „ x 


1=1 


F(x)dx. 


It is clear that 


-'n — k:n 


Hi ~ Hi — 


(F n (x) - F(x))dx - 


F(x)dx. 


In view of Proposition 5 combined with equation (4.9) in Csorgo (1996), we have 
for any x < Z n _ k:n , 

F n (x) -F(x) 


F(x ) 

d (Tty («) - TT™ („)) 

TF ) 


H ’' ( - v) - H{v W'\v)+O r (l/k) 


H 2 (v) 


Integrating the first integral by parts yields 

F n (x) — F (x) _ H i 1} (x) - H [l) _ (x) _ /-gr(i) ^ _ -^( 1 ) ^ 


+ 


F {x) H (x 

r F M - h w w 


H 2 (v) 


dH (v) 


n n (v) H ( v ) dH m (v) + 0p . 
H (v) 


Recall that 

y/n {H n (v) - H (u)) = y/n (v) - H (1) (v )) + y/n (v) - H (0) (v)) , 
which by representations (4.7) and (4.8) becomes 

y/Z (H n (v) - H (v)) = a n (9) - a n (d - H [l) (u)) -a n (l- H (0) (u)) . 
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Furthermore, from the classical central limit theorem, we have (0) — H ^ (0) = 

Op (n -1 / 2 ) . Therefore, we have 

F n (x)-F(x) 1 fi n (x) 1 f x p n (v) 


+ 


F(x) y/nH(x) VnJ 0 h' 2 (v) 


dH (■ v) 


1 f fin T fin ( v ) ,-77(1) ( \ \ /i //-, , D I /~\ 

-/ -zzo- dH (v) + Op (l/re) + Op (1/\M) • 

v B Jo H" (v) 

By letting a n := (k/n) 1 ^ 2 / (hF(ti)) , it is easy to verify that 

Vife (/4i - /ii) 






2=1 


where 


T n i := a n I fid- ’ F(x)dx, 


" Zn - k - n fin(x) 

10 H(x 

^ ^n—k:n 


T n o a„ 


'0 


-'n — fc:n 


Tn 3 • ^r; 


'0 


dH(v) 1 F(x)dx , 
o H (v) J 

/o H (v) 


r^n-k-.n _ 

T„ 4 := a„0 P (y/n/k) / F(x)dx, 

Jo 

T n5 := —a n yfn / F(x)dx and T n6 := 0 P (a„). 


- / n — k:n 


By using the Gaussian approximation (4.10) , we obtain 


(4.11) 


'n—k:n f 

T ri i = a n I =—B n (x) dx 


i o 


tf(x) 


+ op (1) a r , 


2 „-fcT, H (x 


(IB G 1/2 


i/(x) 


-F(x)dx, 


where 


B„ (x) := B n (0) - 5 n (0 - i/ (1) (x)) , for 0 < hf (1) (x) < 0. (4.12) 


Next, we show that the second term of T n i tends to zero in probability, leading to 

r Z n-k:n JPfy) 

T n 1 = a n / =—rB n (x) dx + Op (1). 

do #(u) 
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Let 0 < £ < 1/4 and note that since H = H ^ + H^\ then < H and 


Op (n c ) 


- J n — k:n 


t 7 m ( o 1/2 - c 

H (x 


n — k:n 


- F(x)dx < 0 P (1) n Q a n j 

H{x) Jo 


F(x) 


(H(x)) 


1/2+C 


dx. 


We show that 


n ^a r 


+ (#(*)) 

Indeed, we have 

rZn-k-.n -p( X \ 

>0 


F(x) , _ c f h F(x) 

dx = n ^ a„. I - , . > dx + op (1). 


~ tt / 1/2+C 


' rtit \\ i/2+c 1 


(HM)- 


(H(x)) 


1/2+C 


dx- t - FJ— dx 


0 


(H(x)) 


n — k:n 


F(x) 


, 1/2+C 


dx 


'* (H(x)) 

^max(/i,Z„_fc„) F(x) 


. 1/2+C 


dx. 


J min(h,Z n _ k:n ) (^H (x)) 

By using Potter’s inequalities, given in assertion 5 of Proposition B.1.9 in de Haan and Ferreira 
(2006), we write for e > 0, 


n ^a n 


rmin(h,Z n _ k: , 


F(x) 


iin(h,Z n _ k:n ) (H(x)~) 


1/2+C 


dx < k~ C r X _1/7l + (1/2+C)/7±e l maX B,^n-fc:n/M 

— L' J min(l ,Z n _ k;n /h) 


On the other hand, combining Corollary 2.2.2 with Potter’s inequalities given in 
Proposition B.1.9 (5) in de Haan and Ferreira (2006), yields that Z n _k:n/h —>■ 1 in 
probability. Therefore, the right-hand side of the previous inequality tends to zero, 
as sought. Now, we show that T n] may be rewritten into 


f h F( x ) 

T n i = a n =—(-B n (x) dx + o P (1). 

Jo H(x) 


(4.13) 


Observe that 


f h F(x) r Z n-k-.n jz( x \ 

T n i a n I / cBn (*!') dx T a n I —- (ic) dx T op (1), 

Jo H(x) Jh H(x) 

with the second term in the right-hand side tending to zero in probability. Indeed, 
for fixed 0 < rj, e < 1, we have 


P 


r Z n-k:n f ( V ~\ 

/ ~Trt \ Bn (+) dv 

Jh H(v) 


> T) 


< P 


Z n — 


n—k:n 


h 


-1 


> e + P 


r(l+e)h Jp( v \ 

a n / =-d-B n (v)dv 

Jh H(v) 


>V), 


where, in virtue of the fact that Z n _k:n/h —> 1, the first term tends to zero. It 
remains to show that the second term in the right-hand side is also asymptotically 
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negligible. We have H^' 1 < H, then 


E 


’(i + e ) h f(v) 

B n („ dv 

i H {v) 


™(l+e)/i 


— 


< a r , 


£^Jh { 1) (v)dv 
H 0 ) V 1 ' 


r*(l-\-e)h 


F(v) 


H(v) 


:d,V. 


Changing variables and applying Potter’s inequalities to the regularly varying func¬ 


tion F(x)/\ H(x), yield that, for all large n and £ > 0, we have 


E 


"(l+e)/i 


Hv) 

H(v) 


B n (■ v ) dv 


< a r 


hF{h) 


’l+e 


v -l hl+ l/{2 1 )±i dv 


H(h ) 71 


- 1 +e 


v -i/'n+i/(2n,)±s dVm 


The latter integral is clearly finite and tends to zero as e 4- 0. By similar arguments 
using approximations (4.9) and (4.10) , we also show that 

ph 




— ^-^-dH{y) [• F(x)dx + op (1) 


'o [Jo H (v) 


and 


where 


rh ( px T3 * 


T n :$ — —a r 




I o [Jo H (v) 


f(0) 


dH (v) > F(x)dx + op (1), 


f(o) 


B; (x) := B n (x) - B n 1 - H y ) (x) , for 0 < H y) (x) < 1 - 9. 


(4.14) 


(4.15) 


(4.16) 


Before we examine T n4 , we provide an approximation to T n5 , for which a change of 
variables yields 


r n5 = -Vk 


F(hx) 


dx. 


n — k:n / 


,/h F(h) 

For the purpose of using the second-order condition of regular variation (2.6) of F, 
we write 


T n5 = -VkA^h) 


^n—kin/h' 


F(hx)/F(h)-x~ 1 ^\ , /-f 1 

-—— - dx - Vk x 

A i(h) J Jz n _ k:n /h 


' 17/71 dx. 


(4.17) 

From the inequality (2.3.23) of Theorem 2.3.9 in de Haan and Ferreira (2006), page 
48, we infer that the first integral in (4.17) is equal to 

/■1 


(1 + Op (1)) 


x 


-1/71 f.y.n/71 


— kin/h 


{x Tlhl — l) /71 Tidx, 
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P 

which tends to zero in probability due to the fact that Z n _k :n /h —» 1. Moreover, 
the term \fkA\[h ) has, by assumption, a finite limit. Consequently, the first term 
in the right-hand side of (4.17) is asymptotically negligible. We develop the second 
integral and make a Taylor’s expansion. Knowing, once again, that Z n _ k . n /h —» 1 
ultimately yields that 


r „ 5 = (i + mi)) - iY 

By using result (2.7) of Theorem 2.1 in Brahimi et al. (2015), we get 


T n b — 7a/ ( h ) + op (1). 


(4.18) 


Next, we readily check that the fourth term T n A tends to zero in probability. In- 

pZn-k:n _ _ 

deed, we have / F(x)dx < fi and by assumption \/khF(h ) —> oo. Finally, 

Jo _ 

for the last term T n q we use the second-order regular variation of the tails F and 
G. From Lemma 3 in Hua and Joe (2011), there exist two positive constants c± 
c 2 such that h = (1 + o (1)) Ci (fc/n)~ 7 and F(h) = (1 + o (1)) c 2 (k/n) 1 ^ 1 , thus 
a n = (1 + o (1)) cic 2 (k/n) l ^ 2+1 77,71 . But the indices 71 and y 2 belong to 77, hence 
1/2 + 7 — 7/71 are positive. Therefore, a n —> 0 and T n6 = op (1). The four approxi¬ 
mations (4.13) , (4.14) , (4.15) and (4.18) together with the asymptotic negligibility 
of both T n4 and T n6 give 


Vk (J1 1 - fii) 


hF(h ) 


Ojt. 


l B n (x)- Urn ^ 1 r B > 

-=-7 ( x)dx + a. 


H(x) 


h [ px ■ 

n 1 1 _^ ’ dH(v) F{x)dx 


'0 [Jo 77 (v) 


(4.19) 




a n I < / —-<777 (u) > F(x)dx + 7 W —B* (/i) + o P (1). 

'0 [Jo 77 (n) I V k 


Let us now treat the term \/fc (/i 2 — // 2 ) / (hF(h)) . Consider the following forms of 
/i 2 and / 1 2 : 


/i 2 = hF (/i) 


F (/ix) 


da; and /a 2 = 


Yi 


(H,c) 


1-7) 




Zn—k\nF (Z n — k:n ) 


Fn \Zn—k:n) 

F {Z n _ k:n ) 


OO 
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and decompose yfk (jii — /i 2 ) / (hF (/r)) into the sum of 

7l^’ 6 ) F ( Zn-k-.n ) F n ( Z n _ k , n ) f Z n _ k:r 


Sn 1 ■= Vk 


1 -7l (H ’ C) 


h 


1 , 


5 n 2 := v/fc 


-f 1 (Z n — fc:n) F n i^Z n _ k:n ) I 7i 




7i 


F(/i) F(Z n _ k:n ) ll —7^’ c) l-7i ’ 


S n3 := 


7i F (Z n _ k:n ) j F n ( Z n _ k:n ) ^ 


Sni := 


F n5 := Vk 


l-7i F(/i) 1 F(Z n _, n ) 

7! I F(z n _ k:n ) {Z n _ k .. n ^- lhl 


1 - 71 I F (h) 


7i 


I-71 


Zn—k:n 


h 


-1/71 


h 


1 , 


S n6 := y/k 


1 — 7i ./1 F(ft) 


_ Jr Jr 

For the first term, we have 71 —> 7 ! and Z n _ k:n /h —> 1, which, in view of the regular 
variation of F, implies that F (Z n _ k:n ) = (1 + op (1)) F (h) . Moreover, from (5.29) 
we infer that F n ( Z n _ k:n ) = (1 + op (1)) F ( Z n _ k:n ). It follows that 

S„i = (1 + Or (1)) rr^-Vk _ d , 


1-71 V h 

which, by applying result (2.7) of Theorem 2.1 in Brahimi et al. (2015), is approxi¬ 
mated as follows: 

«a - (1 + «,(1)) (ft). (4.20) 

By using similar arguments, we easily show that 

1 


Sn2 — (1 + Op (1)) 


Vk<^ H ' c) 


7i - 7i h 


(1-71 y 

which, by applying result (2.9) (after a change of variables) of Theorem 2.1 in 
Brahimi et al. (2015), becomes 

(;7f *-'« M*"- ifa. (» * ■^ 


For S n 3 , we have 


S n 3 — (1 + Op (1)) 


7i 


l-7i 


Vk ( F - (Zn ~ k - l) . 
V F (Z n _ k:n ) J 
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Using Proposition 5.1, we have 


Sn3 — (1 + op ( 1 )) 



r h B n (v) 
Jo U 2 (v) 


dH(v ) 




/o H {v) 



(4.22) 


+ (i + 0 r (i)) r 2i - 

1 - 7i 


B„ (h) + op (1). 


For the fourth term, we use the second-order condition (2.6) of F and the fact that 

P 

Z n _ k:n /h ->• 1 to get 


S n A = op i^/kAi (h)j = op (1), asn-> oo. (4.23) 

P 

For S n s, we apply the mean value theorem with the fact Z n _ k:n /h —> 1 to have 

s „ 5 = - (i + Ml)) - iV 

Using, once again, result (2.7) of Theorem 2.1 in Brahimi et al. (2015) yields 


S n 5 — — (1 + op (1)) 
For the last term, we first note that 


7 


t b ; w . 


S n 6 

y/k 


X 




X 


1 - 71 V k 

F ( hx) 


(4.24) 


dx. 


'i ./i F(h) 

Then, by applying the uniform inequality of regularly varying functions (see, e.g., 
Theorem 2.3.9 in de Haan and Ferreira, 2006, page 48) together with the regular 
variation of |7Li|, we show that 

VkA 1 (, h ) 


5, 


77.6 


(7l + Ti ~ 1) (1 - 7l) 

By gathering (4.20) , (4.21) , (4.22) , (4.23) , (4.24) and (4.25) we end up with 

Vk (P2 - /i 2 ) 71 [k 


(4.25) 


hF(h) 

I -71 

V n 1 , 

Jo II' 



71 B 

n(h) 


H (v) - 


, dH y, (v) 

0 H\v) 


r>00 


k j p (l - 7l ) 2 

+ Rnl + Op (1), 


7B; (h) + J -A- 


v (. hv ) dv 


P(1 -7i)" 


(4.26) 


VkA 1 (h) f 1 1 ] 

(! - 7 i) I (1 - P T 1) (1 - 71) (71 + n - l) J ' 


where 
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Finally, by summing up equations (4.19) and (4.26) we obtain 

D »‘ ++mi). 


hF(h) 


2—1 


where 


D n i := a n [ F{v)dv , D n2 := a n [ If dH(v ) l F(x)dx , 


iJ(u) 


'o [do H (v) 


-^n3 • 


B f^W (1) (u) |>F(x)dx, 


'o [do H (v) 


D n A • 


-Dn5 : ~ 


7! [k f h B n (v) 


l-7i V n Wq jr( v ) 


dif(u) 


’b;(u) jT7 (d 


7 i 


k V p(l -7i) 


r B n (h) + 


/o 5"» 

1 ^°° 


^ » , 


p(l-7i) Ji 


v l B* n (hu) dv 


Note that D n2 is of the form —ci n / -0 (x) dip (x) , where <p (x) := / F(u)du and 


,—2 


-0 (x) := / B n (v)/H (v)dH(v). Integrating by parts yields 


r h J 3 _ IF / F(x)dx rh 

D n2 = a n p( v )^-dH(v)-‘l Jh 


Bn(^) ■=, v 


I o iT(u) V n /iF(/i) do H (v) 

Equation (5.1.9) in Theorem B.1.5 (Karamata’s theorem) in de Haan and Ferreira 

poo 

(2006) yields that / F(x)dx/ ( hF(h )) —> 71 / (1 — 71 ). We apply the same tech- 

Jh 

nique to D n3 and get 


D 


n 2 


D 


n3 


D n4: — L n2 + L n 3 + R 


k n2 5 


where R n2 := op(5 n4 ) and 

r f h B n (v) , f h B* n {v) 

Ln2 ■= a n / =2 — 00 «d/(u) and L n3 := -a n / —p (v) dH (v). 

do 5 ( v ) do 5 (u) 

This yields the following new decomposition: 


— ^ + 5 n i + i? ri 2 + op(l), 

with L nl := D n i and L n4 := D n5 . The four L ni are centred Gaussian rv’s whose 
asymptotic second moments are finite, as we will see thereafter. Indeed, L n4 is 
the Gaussian approximation to Hill’s estimator given by result (2.9) of Theorem 
2.1 in Brahimi et al. (2015), hence we have lim^oo E [L^ 4 ] < 00 . For the three 
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others, we literally compute the asymptotic moments of order two. Note that from 
the covariance structure in Csorgo (1996), page 2768, we have the following useful 
formulas: 

E [B n (u) B n (u)] = min (■ u ), H (l) (v)'j - H {1) (■ u) H (1) (■ v ), 

< E [B; (m) B; (v)] = min (if (u), H (u)) -H(u)H (v), (4.27) 

^ E [B n (u) B; (u)] = min (h { 1) (u) , H (1) (u)) - H (1) (u) H ( v ). 


After elementary but very tedious calculations, using these formulas with l’Hopital’s 
rule, we get as n —> oo, 


k 


h ^E[B n (u)B n (u)] ^ 

_ 9 -- dH ( u ) dH ( v) —> p, 

n o Jo H (u)H {v) 





n 



o Jo H ( u)H (v) 





n 7o Jo H (u) H {v) 
By using the results above, we obtain 

E [L n i] —> 

E [7v n2 ] 2 —>■ 


2 7 2 7i 2 


( 7 i - 7 + 77 i) ( 7 i - 27 + 2771) ’ 
2P7 f 


and 


E [L n ^ —> 


(71 - l) 2 (71 - 7 + 771) (71 - 27 + 2771) ’ 

_2p 2 7i_ 

(71 - l) 2 (71 - 7 + 771) (71 - 27 + 2771) 
As a consequence, we conclude that 

\Pk ^= - —%■ J\f (m, a 2 ) , as n —» 00 , 

/iF(/i) V ' 


(4.28) 


where m := lim n ^ 00 .R n i and a 2 := lim^ooE [Y^t=i Lni] 2 ■ The expression of m 
is simple and easily obtainable whilst that of a 2 is very complicated and requires 
extremely laborious computations. However, we readily check that, it is finite. 
Indeed, in addition to the finiteness of the asymptotic second moments E [L ni ] 2 , the 
asymptotic covariances E [L m L nj ] are, in virtue of Cauchy-Schwarz inequality, finite 
as well. Finally, we use the facts that Z n _fc :n /h and F(Z n -km)/F(h) tend to 1 in 
probability to achieve the proof. □ 
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5. Appendix 

In the following basic proposition, we give an asymptotic representation to the 
Kaplan-Meier product limit estimator (in Z n _k:n)- This result will of prime im¬ 
portance in the study of the limiting behaviors of many statistics based on censored 
data exhibiting extreme values. 


Proposition 5.1. Assume that the second-order conditions (2.6) hold. Let k = k n be 
an integer sequence satisfying, in addition to (1.3), VkAj ( h ) = O (1), for j = 1, 2, 
as n —y oo. Then there exists a sequence of Brownian bridges {B n (s); 0 < s < 1} 
such that 


Vk 


(Fn (^Z n —k:n) 

F ( Z n _ k:n ) 






B„, (v) 
H 2 (v) 


dH (v) 



b; (tQ 
H\v) 


dH “> 



+ Op (1) , 


where B n (v) and B* (v) are respectively defined in (4.12) and (4.16) . Consequently, 


Vk ( - l) 4 N (0,p), as n -» oo. (5.29) 

V T \^n-k:n ) J 


Proof. Multiplying (4.11) by \fk yields 


F n fx) — F (x) [k ( e ) 


Vk 1 -^ 


F(x) 


-F7(!) 


H(x) 


k r x a n (#) - Onto- H (1) (v)) _ 

1 V - J -dH(v) 


H 2 (v) 



l~jf y* a n (0) - a n (d - H (l) (v)^j - a n (l - H (0) (u)) 

V Wo H 2 (v) 


(u) 



The Gaussian approximations (4.9) and (4.10) , in x — Z n _ k . ni and the facts that 
a Jk/n and 1 /Vk tend to zero as n —* oo, lead to 


Vk 


Fn \Z n —k:n) F (Z n — k:n ) 


F(Z n _ 


Bn (Zn—k-.n) T 


k\n) 


k r z -^KM dlI iy v)+0rm 

o H (i>) 



B„ (u) 
H 2 (v) 


dH (v) 
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Applying Lemma 5.1 completes the proof. The asymptotic normality property is 
straightforward. For the variance computation, we use, in addition to (4.28) , the 
following results: 


rE[BU«)B„ w W and r" (M) _ Vi 

/o iT (u) do 


H 2 (u) 


similarly obtained as (4.28) . 


□ 


Lemma 5.1. Assume that the second-order conditions of regular variation (2.6) and 
let k := k n be an integer sequence satisfying (1.3). Then 

" Z n-k:nT3 _ 

-%UdH(v) = or( 1). 


I u 

U) \ I — / 

n -'h 


H v) 


r .s fk r Zn - k:n B* n (v) — (i) 

O) \/r I — 2 , dI1 (^) = °p( 1 ) 

H v) 


n -i h 

n 


( Hi ) ^ B n (h)} Op (1) 


H {b; ( z n _ k:n ) - b; (/*)} = o P (i). 

Proof. We begin by proving the first assertion. For fixed 0 < 77 , e < 1, we have 


P 


< P 


- J n — K\n 


B n (v) 


dH (v) 


n Jh 
Zn—k:n 


H 2 (v) 


h 


> 6 + P 


> Tj 

~jf l‘(l+e)h 


dH (n) 
B„ (») =W 


> V I 


H (v) 

It is clear that the first term in the right-hand side tends to zero as n —> 00 . Then, 
it remains to show that the second one goes to zero as well. Indeed, observe that 

iff r( i+ e A 


E 


k f {1+e)h „ , dH (v 

B« [V 


H 2 (v) 


< -"nJ h 


E |B n (v) 


dH (v) 
H 2 (v) 


From the hrst result of (4.27) , we have E B n (u)| < \JH^ (v). Then 


E 


k [ {1+e)h n , s dH (v) 
B n (v) 

n l h H (v) 


< 


1 

n -lh H (v) 


which, in turn, is less than or equal to 

lc 


H [1) (h) _ _ 

n ' yJ \H((l + e)h) H(h ) 




























24 


Since 77 ( h ) = k/n, then this may be rewritten into 


H W (h) 
\ H(h) 


\H((l + e)h) ) 


Since H (1) (h) rs./ pH (/i) and H G then the previous Quantity tends to 

p 1 / 2 ^(1 + e) 1 ^ 7 — 1^ as n —> oo. Being arbitrary, e may be chosen small enough so 
that this limit be zero. By similar arguments, we also show assertion (ii ), therefore 
we omit the details. The last two assertions are shown following the same technique, 
that we use to prove (iv). Notice that, from the definition of B* (v) and the second 
covariance formula in (4.27) , we have 


(B* (v); v>0 } = {B n {H{v))- v > 0} , 


where {B n (s); 0 < s < 1} is a sequence of standard Brownian bridges, ffence 

yf {B; - B = {Bn (H (Z„- t:n )) - B„ (77 (/.))} . 

Let {W n (t); 0 < s < 1} be a sequence of standard Wiener processes such that 
Bn (t ) = W n (f)-fW n (1). Then ^Jk{W n ( Z n _ k , n ) — B* (/i)} equals in distribution 
to 

Jj: ({'w n (77 (Z n _* :n )) - w n (77 (/!)) } - (77 (Z n _ fe:n ) - 77 (/i)} (1)) . 

Using the facts that 77 (h) = k/n and 77 (Z„_fe :n ) /77 (/i) = (1 + op (1)), we get 

(H(Z n _ k:n ) -71(h)) = ^ (IW±4 _ = M l). 

On the other hand \fnjk (W„ (77 (Z n _fc :n )) — W n (77 (h)) j is a sequence of Gauss¬ 
ian rv’s with mean zero and variance 

^(77(Z„_ fc „)-77(fc)) = LU^L-l. 

which tends to zero (in probability), as n —> oo. This achieves the proof. □ 
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